On recent puzzles in the production of heavy quarkonia 
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Recently, several surprising experimental observations in 
the production of heavy quarkonium have been reported. In 
e + e~ annihilation at y/s = 10.6 GeV, Belle Collaboration 
finds that J/ip mesons are predominantly produced in associa- 
tion with an extra cc pair, with a(e + e~ — > J/ipcc)/a(e + e~ — > 
J/i/>X) = 0.591°;^ ± 0.12, and the BaBar collaboration re- 
ports that the produced J/ip's have mostly longitudinal po- 
larization. In pp collisions at the Tevatron, the CDF Col- 
laboration reported an excess of J /if) and ip' mesons at high 
p_L over the perturbative QCD predictions; non-perturbative 
approach of NRQCD can accomodate the magnitude of the 
production cross section but not the observed experimentally 
polarization of quarkonia. In this note we propose possible 
solutions to these puzzles, and devise further experimental 
tests. 



Experimental and theoretical investigations of heavy 
quarkonium have been very important in establishing 
QCD as a true theory of strong interactions. The prop- 
erties of charmonium and bottomonium states have been 
found in qualitative and in some cases quantitative agree- 
ment with theoretical expectations based on QCD. Nev- 
ertheless, since many important non-perturbative prob- 
lems, including confinement and the spontaneous viola- 
tion of chiral symmetry, are still not solved, one may 
expect that a closer look at the properties and the pro- 
duction of quarkonium states will help in elucidating the 
interplay of short and long distances in this theory. This 
is why new experimental results on quarkonium produc- 
tion always attract much interest and attention. 

Recently, Belle Collaboration reported the measure- 
ment of J/if) production in e + e~ annihilation at y/s — 
10.6 GeV [1]. It was found that the cross section of 
J /if) production significantly exceeds theoretical expec- 
tations based on the Color Singlet Model (CSM) [2] and 
non-relativistic QCD (NRQCD) [3,4]. Even more sur- 
prisingly, it was found that most of the observed J/if>'s 
were accompanied by an extra cc pair, with cr(e + e~ — > 
J/ipcc)/a([e + e- -► J/ipX) = 0.59 j^f ± 0.12. This ratio 
exceeds the existing theoretical predictions by about an 



order of magnitude. The problem has recently attracted 
attention of many theorists - see, e.g. [5-12]. 

At first glance, it seems impossible to understand this 
result within the framework of QCD, where the produc- 
tion of an extra cc pair should be suppressed both by 
the large mass of the heavy quarks and by extra powers 
of the strong coupling constant a s . However, this con- 
clusion may be premature. To explain our point, let us 
consider the diagrams of J/ip production in e + e~ anni- 
hilation. 

There are three types of such diagrams. In the dia- 
grams of the first type (Fig. la) J/ip is formed by the 
fusion of cc pair produced by the initial virtual photon. 
This means that cc quarks, which initially were moving 
in the opposite directions, turn around and have almost 
equal and parallel momenta in the final state (up to the 
terms ~ rj/^Ej/^/Mj/^, where rj/$ and Ej/^ are the 
radius and the energy of J/ip). The momentum con- 
servation is ensured by the radiated gluons which may 
produce light (q = u, d, s) qq or charmed cc pairs. The 
diagrams of the second type (Fig. lb) correspond to the 
fragmentation of c (or c) into J/ip. This process requires 
the production of an additional cc pair by emission of at 
least one gluon by initial c or c. Finally, the third type of 
processes arc described by the diagrams in which the ini- 
tial virtual photon creates a pair of light quarks and the 
J / if) is formed from the cc pair produced by an exchange 
of gluons (Fig.lc). Evidently, the diagram of Fig.lc is 
suppressed in comparison to diagrams of Figs. la, lb by a 
factor of a s (m c ) and will be disregarded in what follows. 

On general grounds one may expect that the diagrams 
of Fig. lb dominate at very high energies of the colliding 
e + e~, when the ladder of cc pairs is formed and the pro- 
cess may be described by Regge theory (such process was 
considered by Kaidalov [12]). It is easy to estimate the 
energies starting from which one may expect the domi- 
nance of the diagrams of Fig. lb. Let p c be the momen- 
tum of the c— quark (in the e + e~ cm. system) frag- 
menting into J I if) (see Fig. 2). Then the minimal value 
of the recoil momentum q, corresponding to the forward 
production, is equal to 
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By requiring q to be at least as small as q ~ 0.5 GeV 
(typical for Regge asymptotics), we get p c > 10 GeV, i.e. 
y/s > 20 GeV. In fact one may expect that the energy 
should be y/s > 50 GeV. 

Fig. la dominates at low and intermediate y/s. Our 
main interest here will be concentrated on the qualitative 
and semi-quantitative consideration of its contribution at 
intermediate i/s ~ 10 GeV where the measurements of 
Belle and BaBar were done. 



J/tp. The recoil mass was defined as 
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FIG. 1. The production of J/ip in e + e annihilation. 

We will use the well-established notion of duality be- 
tween the quark-antiquark continuum of pQCD and the 
spectrum of physical bound states, and assume that in 
order to produce the J/ip (or any bound charmonium 
state in general) the integral over the invariant mass y^s^ 
of the cc pair in the diagrams of Fig. 1 should be restricted 
to small values corresponding to the typical velocity of 
charm quark in the J/ip wave function, ~ 0.45. Us- 
ing v = yl — 4m 2 /s c , and the pole mass of the charm 
quark m c ~ 1.3 GeV, one finds that this corresponds to 
the invariant mass y/s^ = 2m c /(l — v^) 1 / 2 ~ Mj/^. 




FIG. 2. The kinematics of fragmentation into J/tp- 

This means that the invariant mass of the hadron sys- 
tem produced in association with J/ip in the e + e~ — ► 
J/ipX process at y/s = 10.6 GeV can be quite large. Belle 
Collaboration [1] measured the distribution of events as 
a function of the mass of the system recoiling against the 
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where E*j^ and Vj^ are the energy and momentum of 
J/ip in the cm. frame. It was found that M recoi i con- 
centrates in the domain M reco u > 5 GeV, with almost 
no events below 2.8 GeV [1]. 

We note that this invariant mass is sufficiently larger 
than the threshold for the production of an additional 
cc pair. Let us now estimate the relative probabilities of 
the associated with J/ip production of heavy and light 
quark-antiquark pairs. 

As the recoiling invariant mass available for the ex- 
tra quark pair production is large, M recoi i > 5 GeV, let 
us, in the first rough approximation, treat c quarks as 
massless. Since three different flavors (q = u, d, s) of 
light quark pairs can be produced, in perturbation the- 
ory one should sum over them incoherently, and we ex- 
pect a{J/ipcc)/a{J/ipX) ~ 1/4 for X = cc + qq. Non- 
perturbative effects, however, can induce correlations in 
the production of quark pairs. In this case, since the qq 
pairs are produced by gluons, the wave function of the 
light quark state Xu g ht has to be an SU(3) singlet: 
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Therefore, light and charm quarks should be produced 
in association with the J/ip with equal probabilities, and 
we find 

R = a(e+e- -» J/ip cc)/a(e + e - -» J/ip X) = 



\{e+c-\J/iP cc)| 2 

\(e+e~\J/ip cc) | 2 + 1 1 (e+e~\J/ip uu + dd + ss) | 2 
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The result (4) corresponds to the assumption of maxi- 
mal coherence in the production of light quarks. Under 
the opposite assumption of completely incoherent pro- 
duction of light flavors we would have R ~ 1/4. Basing 
on the data [1] M recoi i may be estimated as M recoi i w 5.5 
GeV and 4m 2 /M 2 ecoi; w 0.2. Therefore, we estimate m c 
corrections in (4) as O(4m 2 /M 2 ecoi/ ) ^0.2 and our ex- 
pectation for R is 
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which is not far from the Belle result. 

Theoretically the suppression of the spectrum at 
M recoi i < 3 GeV may be understood if we account for 
the fact that the two gluon state recoiling against J/ip 
is a scalar; because of the scale anomaly in this channel, 
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there are strong nonperturbative effects resulting in the 
high mass of the produced gluonic states. 

In addition to the diagram of Fig. la, at sufficiently 
large yfs one should also consider the diagram of charm 
quark fragmentation, see Fig. lb. This diagram always 
produces an extra cc pair in association with the J/tp, 
giving the ratio (4) equal to unity, and the admixture 
of the heavy quark fragmentation will thus increase our 
estimate for the ratio (4). 

Let us now turn to the discussion of J/tp polariza- 
tion. Since the photon in e + e~ annihilation is mostly 
transverse, it has helicity A = ±1. Therefore the cc pro- 
duced by the photon should have opposite helicities of 
A c = +1/2, A s = -1/2 or A c = -1/2, A e = +1/2. Ini- 
tially, heavy quark and antiquark move in the opposite 
directions in the center of mass system of e + e~ anni- 
hilation with the velocities v = yl — 4m;; /s which are 
close to v ^ 1 at y/s = 10.6 GcV. However, to become 
bound in the J/tp (or any other bound state of char- 
monium), at least one of the quarks has to change the 
direction of its momentum by radiating gluons (and ex- 
tra quark-antiquark pair(s)) since the relative velocity of 
heavy quarks in a bound state should be small. Since the 
coupling of quarks to gluons q^^qA^ in QCD conserves 
the helicity of the quark, a change in the direction of its 
momentum should be acoompanied by the spin flip. We 
thus come to the conclusion that in the case of J / tp pro- 
duced at high momentum, the total spin of J/tp should 
have zero projection on its direction of motion, which 
corresponds to the longitudinal polarization. 

This is in agreement with the experimental result of 
the BaBar Collaboration, which states that the angular 
distribution of positively charged lepton decay product 
with respect to the direction of J/tp measured in the CM 
frame is W(0) ~ 1 + a cos 2 6 with a = -0.46 ± 0.21 for 
CM momentum p* < 3.5 GeV, and a = -0.80 ± 0.09 
for CM momentum p* > 3.5 GeV. (In this distribution, 
a = — 1 corresponds to longitudinal polarization, a = +1 
to transverse, and a = indicates no polarization). 

Let us discuss now the case of very high y/s > 50 GeV, 
when the dominant process of J / tp production is the c(c) 
quark fragmentation, see Fig. lb. It is evident that in 
this case any produced J/tp is acompanied by an extra 
cc pair, and therefore the ratio 

cr(e + e- -> J/tPcc)/a(e + e- -> J/tpX) ~ 1. 

As is clear from Fig. lb one of the c quarks has helicity, 
say, A = +1/2; the helicity of the other quark created 
from the vacuum is uniformly distributed. Therefore the 
mean value of a is equal to zero (it is easy to see that the 
cases of a — +1 and a = — 1 have equal probabilities), 
and the produced J/tp is unpolarized. 

These two observations show that the fragmentation 
mechanism cannot describe the data [1], [13] at yfs = 10.6 
GeV, in agreement with our expectations. 



A similar mechanism can be expected to work in J/tp 
and tp' production at high p±. Indeed, the dominant 
mechanism of cc production in this kinematical region 
is the gluon fragmentation [16] ; at high p± 3> 2m c we 
may assume that the gluon is almost on shell and thus 
is transversely polarized. Apart from the overall color 
charge, the process thus becomes analogous to the e + e~ 
annihilation that we considered above, and we expect 
that the same conclusion about the longitudinal polar- 
ization of J/tp and tp' should hold there as well. We thus 
predict that the polarization of J/tp and tp', as well as 
of bottomonium states, will become increasingly longi- 
tudinal with the increase of p± . This observation seems 
to be consistent with the CDF data [14,15], but better 
statistics and a larger range in p± would be needed to es- 
tablish this conclusively. Recently it was suggested that 
the fragmentation of charm quarks can dominate the pro- 
duction of J/tp and tp' at p± [17]; if this is the case, then 
as we argued earlier, no polarization should be present 
at high p t . 
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